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Abstract. After a brief introduction to several variational problems in the 
study of shapes of thin structures, we deal with variational problems on 2- 
f^} dimensional surface in 3-dimensional Euclidian space by using exterior dif- 

ferential forms and the moving frame method. The morphological problems 
| of lipid bilayers and stabilities of cell membranes are also discussed. The 

key point is that the first and the second order variations of the free energy 
determine equilibrium shapes and mechanical stabilities of structures. 
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1. Introduction 



The morphology of thin structures (always represented by a smooth surface M 
in this paper) is an old problem. First, we look back on the history [21]. As 
early as in 1803, Plateau studied a soap film attaching to a metallic ring when the 
ring passed through soap water [25]. By taking the minimum of the free energy 
F = X J M dA, he obtained H = 0, where A and H are the surface tension and 
mean curvature of the soap film, respectively. From 1805 and 1806, Young [37] 
and Laplace [11] studied soap bubbles. By taking the minimum of the free energy 
F = p J dV + X § M dA, they obtained H = p/2X, where p is the osmotic pressure 
(pressure difference between outer and inner sides) of a soap bubble and V is the 
volume enclosed by the bubble. We can only observe spherical bubbles because 
"An embedded surface with constant mean curvature in 3-dimensional (3D) Euclid- 
ian space (E 3 ) must be a spherical surface" [1]. In 1812, Poisson [26] considered 
a solid shell and put forward the free energy F = § M H 2 dA. Its Euler-Lagrange 
equation is V 2 H + 2H(H 2 — K) = [36]. Now the solutions to this equation are 
called Willmore surfaces. In 1973, Helfrich recognized that lipid bilayers could be 
regarded as smectic-A (SmA) liquid crystals (LCs) at room temperature. Based on 
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the elastic theory of liquid crystals [4], he proposed the curvature energy per unit 
area of the bilayer 



'lb 



(k c /2){2H + c ) 2 + kK, 



(1) 



where k c and k are elastic constants. K and cq are Gaussian curvature and spon- 
taneous curvature of the lipid bilayer, respectively. Starting with Helfrich's cur- 
vature energy (1), the morphology of lipid vesicles has been deeply understood 
[15, 24, 28]. Especially, the free energy is expressed as F = p J dV + § M {\ + 
£ib)dA for lipid vesicles, and the corresponding Euler-Lagrange equation is [22]: 



p - 2XH + k c V 2 (2H) + k c (2H + c )(2H 2 - cqH - 2K) = 0. 



(2) 



For an open lipid bilayer with a free edge C, the free energy is expressed as 
F = f M {\ + Sib)dA + 7 § c ds, where 7 is the line tension of the edge. The 
corresponding Euler-Lagrange equations are as follows [3, 31]: 



k c {2H + c )(2H 2 - c H - 2K) - 2XH + k c V 2 {2H) = 

[k c (2H + c ) + kk n ]\ c = 



-2k c — + jk n + k-^ 

OG2 as 



-j{2H + c ) 2 + kK + \ + ~fk g 



c 



c 



(3) 
(4) 

(5) 
(6) 



where k n , k g , and r g are normal curvature, geodesic curvature, and geodesic tor- 
sion of the boundary curve. The unit vector e2 (see also Fig. 1) is perpendicular to 
tangent vector of edge C and normal vector of surface M. Above four equations 
are called the shape equation and boundary conditions of open lipid bilayers. The 
boundary conditions are available for open lipid bilayers with more than one edge 
because the edge in our derivation [31] is a general one. 

Secondly, we turn to the puzzle about the formation of focal conic structures in 
SmA LCs. As we imagine, the configuration of minimum energy in SmA LCs 
is a flat layer structure. But Dupin cyclides are usually formed when LCs cool 
from isotropic phase to SmA phase in the experiment [7]. Why the cyclides are 
preferred to other geometrical structures under the preservation of the interlayer 
spacing [2]? This phenomenon can be understood by the concept that the Gibbs 
free energy difference between isotropic and SmA phases must be balanced by 
the curvature elastic energy of SmA layers [17]. The total free energy includes 
curvature energy, volume energy and surface energy. It is expressed formally as 
F = § S(H, K, t)dA, where t is the thickness of the focal conic domain; H and 
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K are mean curvature and Gaussian curvature of the inmost layer surface, respec- 
tively. The Euler-Lagrange equations corresponding to the free energy are as fol- 
lows [19]: 

£{d£/dt)dA = (7) 

(V 2 /2 + 2H 2 - K)dS/dH + (V • V + 2KH)dS/dK - 2H£ = 0. (8) 

Solving both equations can give good explanation to focal conic domains [19]. The 
new operator V can be found in in the appendix of Ref. [32]. 

Thirdly, let us see carbon related structures. There are three typical structures com- 
posed of carbon atoms: Buckyball (Ceo), single-walled carbon nanotube (SWNT), 
and carbon torus. In the continuum limit, we derive the curvature energy [23, 30] 
of single graphitic layer E = J \k c (2H) 2 + kK dA from the lattice model [13], 
where k c and k are elastic constants. The total free energy of a graphite layer is 
F = J ^k c (2H) 2 + kK dA + X J dA, where A is the surface energy per unit 
area for graphite. Please note that the surface energy per unit area for solid struc- 
tures is not as a constant quantity as the surface tension for fluid membranes. The 
Euler-Lagrange equation corresponding to the free energy is V 2 H + 2H(H 2 — 
K) — XH/k c = 0. Ceo and carbon torus can be understood with A = 0, while 
SWNT satisfies R 2 = k c /2X, where R is its radius. 

The rest of this paper is organized as follows. In Sec. 2, we show how to derive 
the Euler-Lagrange equation from the free energy functional by using exterior dif- 
ferential forms. The method is developed in Refs. [31] and [32], which might be 
equivalent to the work by Griffiths [9] in essence. But it is more convenient to 
apply in variational problems on 2D surface in E 3 . In Sec. 3, we give several an- 
alytic solutions to the shape equation of lipid vesicles, and to the shape equation 
and boundary conditions of open lipid bilayers. In Sec. 4, we discuss the elasticity 
and stability of cell membranes. A brief summary is given in the last section. 



2. Variational problems on 2D surface 

Many variational problems are shown in Introduction. Here we deal with them by 
using exterior differential forms. 

Let us consider a surface M with an edge C as shown in Fig. 1 . At every point P 
in the surface, we can choose a right-handed, orthonormal frame {ei, e2, 63} with 
e3 being the normal vector. For a point in curve C, e± is the tangent vector of C 
such that e2 points to the side that the surface is located in. 
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The differential of the frame is expressed as 

dr = L0%ei + u; 2 e2 (9) 
dej = ujijej, (i = 1, 2, 3) (10) 

where U\,ui2, oJij = —ojji = 1, 2, 3) are 1 -forms. Please note that the repeated 
subindex represents the summation from 1 to 3, unless otherwise specified in this 
paper. 



e 3 




Figure 1. A surface M with an edge C. { ei , e2 , e3 } is an orthonormal 
frame. 



The structure equations of the surface are as follows: 

d(Ji = UJ\2 A U>2 (11) 
du>2 = 0^21 A LOl (12) 

o>i3 = auji + bu>2, 0J23 = bwi + cw 2 (13) 
du>ij = Uik A u>kj (i,j = 1,2,3), (14) 

where a, b, c are related to the curvatures with 2H = a + c and K = ac — b 2 . 
The variation of the frame is denoted by 

5r = Sir + 5 2 r + 5 3 r (15) 
8 i r = n i e i (a = 1,2, 3) (16) 
Siei = Qujej, = 1,2,3), (17) 

with = —Qiji I = 1,2, 3). In equation (16), the repeated subindex does 
not represent summation. It is easy to prove that the operator 5i(l = 1, 2, 3) has the 
similar properties with the partial differential while the operator S has the similar 
properties with the total differential operator [32]. 
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Using d(5/r = <5/dr and d^e^ = <5/dej, we obtain variational equations of the frame 
as follows [32] : 

SiUi = dS7i — u^Slm (18) 

5lU)2 = £l\U\2 - ^1^112 (19) 

Oii3 = aQ 1 , Oi23 = &^i; (20) 

62^1 = ^2^21 — ^2^221 (21) 

£2^2 = - ^1^212 (22) 

O213 = bQ 2 , ^223 = C^2! (23) 

(5 3 o;i = O3W31 - ^2^321 (24) 

5 3 u;2 = ^3^32 - ^1^312 (25) 

d^ 3 = O313CJ1 + 1^323^2 ; (26) 



SicJij = 


dQuj + QukLOkj — UikQikj- 


(27) 


Using them, we can prove that 






5i(£ui A 0J2) '■ 


= d(£cj2^i) 


(28) 


5 2 {£uji Aw 2 ) : 


= — d(£ujiQ2) 


(29) 


S 3 ((jJi Auj 2 ) ■■ 


= -2HVL3UJ1 A UJ2 


(30) 


Ss(2H)coi A L02 ■ 


= 2(2H 2 - K)Sl 3 ui A uj 2 + d * df} 3 


(31) 


53KLO1 A UJ2 : 


= 2KHVL3UJ1 AW2 + d*d^ 3 , 


(32) 



where £ is a function of 2i? and * is Hodge star operator [35] satisfying 
the following properties: (i) */ = fu\ A UJ2 for scalar function /; (ii) *uj\ = 
uj 2 , *ll>2 = — oj\. * and d are new operators defined in Ref. [32] that satisfy: 
(i) If df = /iwi +J2U2, then df = fiui 3 _ + /2W23; (ii) *^13 = ^23, *f 23 = 
-wis; (iii) J M (fd*dh - hd*df) = J dM (f*dh - h*df), J M (fdidh - hdidf) = 
Jg M (f*dh — h*df) for any smooth functions / and h on M. 1 

Now, we consider the variational problem on closed surface. In this case, the 
general functional is expressed as: 

F= I £(2H,K)dA+p I dV. (33) 

JM JV 



^These two expressions are very similar to the second Green identity. Their proof can be found in 
Lemma 2.1 of Ref. [32]. 
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Using Stokes theorem and the variational equations of the frame, we can prove that 



dV 



v 



dV = 



v 



. I dV = I : 

JV JM 



n 3 dA. 



(34) 
(35) 



Combining them with equations (28)-(32), we have 8\F = 5 2 F = 0, and the 
Euler-Lagrange equation corresponding to functional (33): 



(V 2 + AH 2 - 2K 







2H 



£+p = 0. (36) 



d{2H) 

Next, we consider the variational problem on open surface with an edge C. In this 
case, the general functional is expressed as: 



F=[ £{2H,K)dA + I T(k n ,k q )ds. 

JM JC 



Similarly, we derive its Euler-Lagrange equations as 

d£ B£ 
(V 2 + AH 2 - 2K)^^ + (V • V + 2KH)— -2H£ = 



(37) 



(38) 



e 2 • V 



d£ 



d(2H) 
d ( dV 



~ ( d£ \ d ( d£\ d 2 ( dT 

+ e2 ' v [MJ ~ Ts{ T9 dKj + dT 2 {dk- n 



or 



d 



dT 



dT 



+ Tn-r- \sr- + -7- r„— - T - —k„ k n 



ds \ dk, 



ds 



dk n 



dk g 9 



c 



d£ 



d(2H) 

<P_ ( dT 

ds 2 \ dk„ 



d£_ 9T 

n 8K + dk g n 



+ K 



dT 

dk n 



d f dT 



ds \dk 



d_ 

ds 



dT 

'dkZ. 



dk n 

^k^ 
dk„ 9 . 



= 



(39) 
(40) 



c 



dT 

+ 2(k n -H)k g — 



0. 



(41) 



c 



In special cases, above equations (36), (38)-(41) are degenerated into different 
equations mentioned in Introduction. For example, if taking £ = En + A and 
T = 7, These equations are degenerated into equations (3)-(6), respectively. 

3. Morphology of lipid bilayers 

There are three typical solutions to shape equation (2): Sphere, torus, and bicon- 
cave vesicle. 

First, a spherical vesicle with radius R satisfies equation (2) if pR 2 + 2\R — 
k c co(2 — cqR) = is valid. Secondly, Ou-Yang used equation (2) to predict that a 
torus with the radii of two generating circles r and R satisfying R/r = \[2 should 
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be observed in lipid systems [20]. This striking prediction has been confirmed ex- 
perimentally by three groups [14, 16, 27]. Thirdly, the first exact axisymmetric 
solution with biconcave shape, as shown in Fig. 2, was found under the condition 
of p = A = as [18]: 

rp 

z = z + tamjj(p')dp' (42) 
Jo 

s'm^(p) = -c pln(p/p B ),c > 0, (43) 

where z(p) is the contour of the cross-section, z axis is the rotational axis, and 
ip(p) the tangent angle of the contour at distance p. This solution can explain the 
classic physiological puzzle [8]: Why the red blood cells of humans are always in 
biconcave shape? 



z 




ft p 



Figure 2. Schematic figure of the contour of biconcave vesicle in the 
first quadrant. 

To the shape equation (3) and boundary conditions (4)-(6) of open lipid bilayers, 
we can find two analytical solutions [31]: One is the central part of a torus and 
another is a cup-like membrane shown in Fig. 3. Numerical method and solutions 
to these equations can be find in Ref. [33]. 




Figure 3. The centtal part of a torus (left) and a cup-like membrane 
(right). 
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4. Elasticity and stability of cell membranes 

A cell membrane is simplified as lipid bilayer plus membrane skeleton. The skele- 
ton is a cross-linking protein network and joint to the bilayer at some points. We 
know that the cross-linking polymer structure also exists in rubber at molecular 
levels. Thus we can transplant the theory of rubber elasticity [29] to describe the 
membrane skeleton. 

Based on Helfrich's theory and physics of rubber elasticity, the free energy of a 
closed cell membrane can be expressed as [32] : 



with E H = (k c /2)(2H + c ) 2 + A and 8 d = (k d /2)[{2J) 2 - Q], where 2J = 
£11 + £22, Q = £\\£22 — e?2- Here = 1, 2) represents the in-plane strain 

of the membrane, k d is the elastic constant representing the entropic elasticity [5] 
of membrane skeleton. 

Using the method in Sec. 2, we obtain the shape equation and the in-plane strain 
equations of the cell membrane as [32] : 



fc d [-d(2J) A uj 2 ~ -(enduj2 - £12^1) + -d(ei 2 ^i + £22^2)] = ( 46 ) 
k d [d(2J) Awi- ^(ei 2 du;2 - £22^1) - ^d(£nwi + £12^2)] = 0. (47) 

An obvious solution is the spherical cell membrane with homogenous strains: 
£11 = £22 = £ (a constant) and £12 = 0. The radius R of the sphere must sat- 
isfy 



Now we will show the biological function of membrane skeleton by discussing 
the mechanical stability of a spherical cell membrane. Using Hodge decomposed 
theorem [35], Cli and $^2 can be expressed as f2io>i + ^2^2 = dQ, + *d\ by two 
scalar functions VL and x- Through complex calculations, we obtain the second 




(44) 




(45) 



pR 2 + (2A + 3k d e)R + k c c (c R - 2) = 0. 



(48) 
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order variation of the free energy for spherical membrane b 2 T = G\ + G2, where 

Gi= <f nl{3k d /R 2 + (2k c c /R 3 ) +p/R}dA 

JM 

+ / n 3 V 2 n 3 {k c c /R + 2k c /R 2 + P R/2}dA 

JM 

+ / k c (v 2 n 3 ) 2 dA + ^ / n 3 v 2 ndA 

JM R JM 

+ kA (v 2 n) 2 dA + J^ I W 2 ndA, (49) 
JM v ' 2R Z J M 

and G2 = (k d /4:) § M V 2 x(V 2 X + 2x/ J R 2 )dA Because G2 is positive definite, 
we merely need to discuss G\. and SI in the expression of G\ are arbitrary 
functions defined in a sphere and can be expanded by spherical harmonic functions 
[34]: fl 3 = E£ TZ=-i ai m Yi m (9, 0) and Q = ££ YZ=-l bi m Y lm {6, (f>) with 
a lm = (-l) m oj ,-m and b\ m = {-l) m h _ m . Considering (48), we write Gi in a 
quadratic form: 

00 1 

d = J2 E 2 l«/™! 2 { 3 ^ + P(* + 1) - 2] + l^/i? 2 - fc c c /i? - 

Z=0 m=0 



r> 

1=0 m=0 

00 i 



+ E E % [ 2l2 ( l + l ? " Kl + 1)] \bi m \ 2 - (50) 



Z=0 m=0 



It is easy to prove that, if p < Pl = m **f_ 1]R + 2k ^]' c ^ (I = 2, 3, ■ ■ ■ ), 
then G\ is positive definite. Thus we must take the minimum of pi to obtain the 
critical pressure 2 : 

1 vp; + ^ (1 _ 2cojR); {3kdR 2 > 121kc) 

Taking typical data of cell membrane, k c ~ 20&bT [6], fc^ ~ 6 x W~ 4 kBT/nm 2 
[12], ~ 4nm, i? ~ 1/xm, cqR ~ 1, we havep c ~ 2 Pa from equation (51). If not 
considering membrane skeleton, that is kd = 0, we obtain p c ~ 0.2 Pa. Therefore, 
membrane skeleton enhances the mechanical stability of cell membranes, at least 
for spherical shape. 



9 

In Ref. [32], we ignore the effect of in-plane modes fii and Q.2 on the critical pressure and obtain 
the invalid value. 
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5. Summary 

In above discussion, we introduce several problems in the elasticity of biomem- 
branes, smectic-A liquid crystal, and carbon related structures. We deal with these 
variational problems on 2D surface by using exterior differential forms. Elasticity 
and stability of lipid bilayers and cell membranes are calculated and compared with 
each other. It is shown that membrane skeleton enhances the mechanical stability 
of cell membranes. 
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